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Employing a spin-parameterized Kleldysh-Usadel technique for the diffusive regime, we demonstrate that 
even in the low proximity limit, considerable long-ranged triplet supercurrents can be effectively generated 
by spin-singlet supercurrents flowing parallel to the interfaces of uniform double ferromagnet interlayers with 
noncollinear exchange fields “ independent ” of actual junction geometry. The triplet supercurrents are found to 
be most pronounced when the thicknesses of the ferromagnet strips are unequal. To experimentally verify this 
generic phenomenon, we propose an accessible and well-controlled structure that can fully isolate the long-range 
triplet effects. 

PACS numbers: 74.50.+r, 74.25.Ha, 74.78.Na, 74.50.+r, 74.45,+c 


Spin carriers play a crucial role in spintronics devices where 
the spin of the carriers are used to transport or store infor¬ 
mation [1], The spin current in diffusive ferromagnets due 
to spin-polarized charged carriers has a short-range propaga¬ 
tion and limited spin coherence, thus presenting a detriment 
to the functionality of most spintronics devices [2], This ob¬ 
stacle, however, can be resolved within the cryogenic realm 
by utilizing ferromagnetic (/■’(/superconductor (5) hybrids to 
generate spin polarized superconducting correlations that can 
propagate over several hundred nanometers with limited de¬ 
cay [3-6]. The induced spin-triplet odd-frequency supercon¬ 
ducting correlations were theoretically predicted [6] for sys¬ 
tems comprised of an s-wave superconductor in proximity to a 
ferromagnet. The signatures of such long-range proximity su¬ 
perconducting spin correlations were found shortly thereafter 
in experiments [7], Numerous works since then, both experi¬ 
mentally and theoretically, have been devoted to the study and 
generation of superconducting correlations with net spin [8- 
12 ]. 

It has been shown that it is possible to study long-range su¬ 
percurrents in ballistic bilayer [13] and diffusive trilayer [8] 
magnetic Josephson junctions where the supercurrent flows 
transversely relative to the noncollinear F/F interfaces. The 
situation is unfavorable however in the low proximity limit 
[14] where only faint signatures of the triplet supercurrents 
with spin projection m=± 1 on the quantization axis exist, 
even in uniformly magnetized trilayers [8, 14]. The low prox¬ 
imity limit is realized for low transparency S/F interfaces, 
strongly magnetized layers, thick F barriers, and temperatures 
near the superconducting critical temperature, which prevails 
in many experiments [4, 6], Thus, it is of considerable inter¬ 
est to determine the most simple and generic situations that 
permit generation of long-range triplet supercurrents in the 
low proximity limit. Moreover, the findings in this limit may 
pave the way for more pronounced triplet generation in sim¬ 
ilar structures within the ballistic and full proximity limit of 
diffusive regime [4, 6]. 

In this Letter, we demonstrate that by considering singlet 


supercurrent flow “parallel” rather than “transverse” to the 
F/F interfaces, substantial long-range triplet supercurrents 
can be induced regardless of actual configuration and ge¬ 
ometry [9, 12], We make use of the two-dimensional qua- 
siclassical Keldysh-Usadel approach in configuration-space, 
which is a general formalism appropriate for inhomogenous 
diffusive F/S heterostructures containing generic magneti¬ 
zation and external magnetic field profiles [15, 16]. This 
two-dimensional approach is capable of explorations previ¬ 
ously unattainable in strictly one-dimensional systems. A 
spin-parametrization technique is also incorporated within 
our method to pinpoint the behavior of the odd- and even- 
frequency components of the total supercurrent. To demon¬ 
strate the generality of our findings, we consider a variety of 
structures [9, 12] and specifically establish clear insights into 
the spatial profiles of singlet/triplet supercurrents in three dif¬ 
ferent magnetic structures that host parallel supercurrent flow: 
A noncollinear double /j /F 2 interlayer that is either verti¬ 
cally [Fig. 1(a)] or horizontally [Fig. 1(b)] stacked between 
two S electrodes, or one that is horizontally sandwiched be¬ 
tween a superconducting electrode and a finite sized normal 
metal [Fig. 1(c)]. In order to have supercurrent flow parallel 
to the F 1 /F 2 interfaces in the two latter structures, we intro¬ 
duce an externally applied magnetic field normal to the junc¬ 
tion plane. The structures reveal that short-range spin-singlet 
supercurrents flowing parallel to the F/F interfaces contain¬ 
ing noncollinear magnetizations and unequal thicknesses are 
effectively converted into long-ranged spin-triplet supercur¬ 
rents. The simple structure shown in Fig. 1(c), permits exper¬ 
imental isolation of the long-range triplets and confirmation 
of this generic effect. 

A general three-dimensional quasiclassical formalism for 
diffusive magnetic heterostructures subject to an external 
magnetic field is described by the following Usadel equations 
[15-17]: 

D[d, G[d, G]] + i[e + diag[ft • f ,(h- f) 2 ],G] = 0, (1) 

in which the Pauli matrices constitute the components of 
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FIG. 1. Schematic of the proposed experimental setups of the basic F/S hybrid structures supporting supercurrent parallel to the F 1 /F 2 
interfaces. The ferromagnetic layers, F| and Fo have uniform magnetizations and N represents a nonmagnetic normal metal with a finite 
thickness cV. Superconductivity is induced via s-wave superconducting terminals of infinite extent (labeled Sj. We denote the exchange field 
of each F layer by hi=(h /, hV, /if )=/io(sin /3i cos ai, sin /3i sin a;, cos /3i), in which /3i and cti are spherical angles for i= 1,2. The junctions 
are located in the xy plane and the S/F interfaces are along the y-axis. In panel (a) the interface of the double ferromagnet region is oriented 
perpendicular to the S/F interfaces creating a two-dimensional hybrid. The widths of the rectangular F strips Wp 1 , and Wf 2 are generally 
different, while their lengths, dp, are equal. In (b) and (c) all interfaces are parallel to the y axis. The widths of the rectangular F strips, Wf, 
are the same while the length of each magnetic layer, dp 1 and dp 2 , are not necessarily equal. To generate a supercurrent flowing along the 
Fi / F 2 interface, an external magnetic field H is directed along the z-axis, normal to the junction plane. 


r. We denote the diffusion constant of the medium by I) 
and G, h are functions of coordinates r=(x,y,z) [15, 16]. 
Here, the exchange field describing the ferromagnetic re¬ 
gion, h(r)=(h x (r),h v (r),h z (r)), can take arbitrary direc¬ 
tions. We have defined a 4x4 partial derivative matrix, 
d, as <9= VI— ieA(r)p 3 , in which A is the vector poten¬ 
tial of the applied magnetic field, H. The energy, e, of the 
quasiparticles is measured from the Fermi surface ep. The 
resultant sixteen coupled complex partial differential equa¬ 
tions should be supplemented by the appropriate boundary 
conditions to properly capture the electronic and transport 
characteristics of F/S hybrid structures. We employ the 
Kupriyanov-Lukichev boundary conditions at the S/F inter¬ 
faces [18] and the induced proximity correlations are tuned 
by the barrier resistance £: £( GdG ) • u=[GbcS) G], in which 
n is a unit vector normal to the interfaces and Gbcs is the su¬ 
perconducting bulk solution. Under equilibrium conditions, 
the vector current density is expressed as an integration of 
the Keldysh block, J(r)=Jo / (feTr {/?3 (G[<9, G]) A }, where 
Jo=iVoeD/4, Nq is the density of states at the Fermi surface, 
and e is the electron charge. To gain insight into the triplet 
supercurrent components, we use the spin-parametrization 
scheme for the Green’s function [ 6 ], Therefore, the anoma¬ 
lous component of the Green’s function takes the follow¬ 
ing form in terms of the even (S) and odd (T) frequency 
parts: F(r, e)=i [S(r, e)+T(r, e).r] t v , where T (r,e) = 
(T x ,T y ,T z ), and T x ,T y have m=±l projections along the 
quantization axis, while T z has m =0 [ 6 ]. If we now substi¬ 
tute this decomposition into the Usadel equation, Eq. (1), we 
are able to separate the contributions S, T x , T y , and T z to the 
charge supercurrent. To study the precise behavior of individ¬ 
ual components constituting the Josephson current, we intro¬ 
duce the following decomposition scheme for the current den¬ 
sity: J x ,y=Jx%+J™y+J^ y +Jx%, where Jf y encompasses 
the § terms, and J!p y the T 7 terms (designating 7 = 0 :, y,z). 
Our formalism allows for arbitrary magnetization orientation 
in each ferromagnetic wire. However to most effectively re¬ 
alize the long-range behavior of the decomposed charge su¬ 


percurrent, we take the magnetizations in the F\ and F 2 wires 
to be orthogonal [ 6 ], with exchange fields /ii=( 0 , h v , 0 ), and 
/i 2 = ( 0 , 0, h z ) or equivalently; /3i=ai=n/2, and / 32 = 0 . We 
focus here on uniform F] j F 2 bilayers (the least and simplest 
layered structure) with Fi and F 2 having unequal thicknesses 
as this results in more pronounced singlet-triplet conversion 
[14]. The magnitude of the exchange field in each ferromagnet 
is typically set at |V|=10Ao, the temperature T corresponds 
to T=0.05T C , and £=4. The energies are normalized by the 
superconducting gap at zero temperature, Ao, and the lengths 
by the superconducting coherence length, £, 5 . 

To have absolute comparisons, we constrain the rectan¬ 
gular F strips to have equal total lengths and widths, i.e., 
dp=dpi~\-dp 2 = 2 . 0 ^s and Wp=Wpi-\-]Vp 2 = 2 . 0 £si respec¬ 
tively (see Fig. 1). The S electrodes are given a macro¬ 
scopic phase difference of \{p\=p/2, which is the phase dif¬ 
ference corresponding to the critical supercurrent. In Fig. 2 
we show results corresponding to the setup in Fig. 1(a) with 
IT Vi =0.3 £s, 14'f2=1T£,s, and no applied magnetic field. 
First, the spatial map of the total maximum charge supercur¬ 
rent density shows that the current in the vicinity of the F-\ / F 2 
junction has a nonzero //-component [19] (as mentioned ear¬ 
lier, y/£s= 0.3 coincides with the width of the F\ nano-wire, 
I'l'Vi =0.8V). This result can be compared with Fig. 3(a) for 
an S/F-[ /F’ 2 / 6 'junction with transverse supercurrent flow rel¬ 
ative to the F 1 /F 2 interface (e.g.. Fig. 1(b) with H =0). Fig¬ 
ure 2 reveals that the induced //-component is present over the 
entire junction width as exhibited by the curved quasiparticle 
current trajectories. The spatial behavior of the triplet and sin¬ 
glet contributions to the supercurrent is also shown as a func¬ 
tion of y (along the junction width) at four x locations along 
the junction length. The top and bottom set of panels corre¬ 
sponds to the current in the x and y directions respectively. 
The singlet and triplets components to the supercurrent both 
play important roles in the supercurrent: The amplitude of J y ° 
is comparable to \J y y +J y z \■ These components of J y , and 
hence J y itself, vanish at //=(), 2.0V, corresponding to the 
vacuum borders. The vector plot of the supercurrent density 
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FIG. 2. Spatial map of the critical cur¬ 
rent density J(x,y) and its even- and 
odd-frequency components. The con¬ 
figuration of the S/F 1 /F 2 /S Josephson 
junction considered here is shown in Fig. 
1(a). The junction length dp is equal to 
2 .0£s and the width of each ferromag¬ 
netic strip labeled F\ and F 2 are unequal, 
i.e., 1Tfi=0.3(;s and Wf2 =1-7 £s. Ver¬ 
tical lines in the panels separate the 
two ferromagnetic regions along the y 
direction. 



in Fig. 2 reveals that the current flow in the middle of the junc¬ 
tion (x=£s) is directed entirely along the x-direction through¬ 
out the junction width. The ^-component of the singlet con¬ 
tribution to the current, becomes vanishingly small when 
approaching this central region. On the contrary, the triplet 
contributions J* v and ./® z are maximal there, demonstrating 
optimal singlet-triplet conversion. Another important aspect 
of this type of junction is seen in the behavior of and 
.J yz as a function of y: these triplet components generated in 
one F region penetrate deeply into the adjacent F segment. 
Therefore, two important phenomena arise: First, a singlet 
supercurrent flowing parallel to the interface of the f j / F- 2 bi¬ 
layers is converted into a triplet supercurrent of both spin pro¬ 
jections. Second, there is a deep penetration of odd-frequency 
triplet supercurrents laterally (along y) into the F wires pos¬ 
sessing orthogonal magnetizations relative to the spin projec¬ 
tion of the triplet supercurrents. The signatures of long-ranged 
proximity effects on the critical supercurrents flowing across 
similar structures as Fig. 1(a) are studied in Ref. 12. The 
planar structures with domain wall patterns are also investi¬ 
gated in Ref. 9 where the textured magnetizations generate 
long-ranged supercurrents regardless of either ‘transverse’ or 
‘parallel’ transports [6], We note that the results of Refs. 12 
and 9 also fully confirm the generic scenario introduced here. 

Next we illustrate an S/F\/F 2 /S junction depicted in 
Fig. 1(b). In Fig. 3(a), the external magnetic flux is absent, 
while in (b) a magnetic flux < T>=3 < F»o is applied to the sys¬ 
tem. The geometric dimensions correspond to dpi=0.3^s, 
dF2=l-7£s, and the junction width, WV=2.0£s. In the ab¬ 
sence of an external magnetic field, the charge supercurrent 
has no component along the y direction [15, 16, 20], and is 
constant along x (the supercurrent flows transverse to the F/F 
interface [14]). The relevant charge supercurrent density com¬ 
ponents, J®°, J* y , and J® 2 , thus only vary spatially in the x 
direction and their behaviors are substantially different from 
their counterparts where the supercurrent flows parallel to the 
Fi/F 2 interface [Fig. 1(a) setup] shown in Fig. 2. It is evi¬ 
dent that J* y disappears in F-> while J ® 2 is zero inside the b\ 
segment since we have the exchange fields fti=(0, h v , 0), and 
ft 2 =(0, 0, ft 2 ) acting effectively as triplet spin filters. In other 
words, J sy is generated in F\, becomes localized there and 
then at the F\ / F 2 interface, converts to the short-ranged J sz 


in F '2 (or vice versa). The triplet supercurrent generated in 
one F is closely linked to the local magnetization texture and 
therefore does not penetrate into the other F whose magneti¬ 
zation is orthogonal to the spin direction [13, 14, 19]. Turning 
now to Fig. 3(b), a magnetic field is applied, corresponding 
to a magnetic flux of <1>=.3'1>(,. This induces proximity vor¬ 
tices that result in a nonuniform supercurrent response vary¬ 
ing in both the x and y directions [15, 16, 20], The corre¬ 
sponding supercurrent constituents, Jy°(x,y), Jy V (x,y), and 
Jy Z {x,y), are thus also nonzero along y. The amplitude of 
J%. v {x,y), is smaller than the other components due in part 
to the small F\ width, dpi =0.3£s, consistent with Fig. 3(a). 




FIG. 3. Decomposed components of critical supercurrent density 
in an S/F 1 /F 2 /S junction depicted in Fig. 1(b). The two ferro¬ 
magnetic strips have unequal lengths: dFi=0.3£s, and <(f2=1-7£s, 
while their widths are equal: Wfi=Wf 2 =2.0£s- In part (a), there 
is no applied magnetic field, whereas in part (b) the external mag¬ 
netic field corresponds to a flux of 4>=34>o. The critical current 
components in both cases, (a) and (b), are plotted as a function 
of ®-position at four different locations along the junction width: 

y=0.25£ s , 0.50&, 0.75&, 1.0£s. 
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FIG. 4. Spatial behavior of decomposed diamagnetic supercur¬ 
rent as a function of lateral position inside the S/F 1 /F 2 /N junction 
depicted in Fig. 1(c). Three magnetization orientations of F\ are 
shown: fh=0, 0.3-7T, and 0.5-7T (ai=7r/2), while F 2 remains magne¬ 
tized along 2 (/?2=0). The thickness of the left F\ layer is fixed at 
dpi=0.1(;s, while the thickness of the right F 2 layer and N metal 
are cZf2=L9£s and djv=2.5£s. 


thick N region with minimal decay [14] while the singlet and 
short-range triplet supercurrent components nearly vanish by 
the time they reach the N layer. Therefore, this pure long- 
range triplet current, which is odd in frequency, can be ex¬ 
perimentally isolated in very simple S/F%/F^/N structures. 
Considering todays’ technological advancements, the experi¬ 
mental investigation of the addressed phenomena in this Letter 
are readily accessible [21], 
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ported in part by ONR and by a grant of supercomputer re¬ 
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As seen, J^ z (x, y), is drastically suppressed at the F 1 /F 2 in¬ 
terface, similar to the case in Fig. 3(a). Examining the sin¬ 
glet and triplet components to the supercurrent flowing in the 
y-direction, we see that their contribution has increased dra¬ 
matically compared to propagation in the x direction (nearly 
an order of magnitude or more). The overall spatial behav¬ 
ior of the even-frequency triplet components are different: 
The component with orthogonal spin projection to the mag¬ 
netization direction Jy V [x , y) penetrates now extensively into 
the ferromagnetic regions compared with its ^-component 
counterparts and the component with parallel spin projection, 
Jy Z (x,y). We note that the magnetization direction in F2 is 
orthogonal to the spin polarization of the J sy supercurrent 
component. By considering these findings together with those 
of a S/F1/F2/S configuration where the Fi/F 2 interface is 
orthogonal to the S/F interface in Fig. 2, and the results of 
Refs. 12 and 9, one concludes that: A singlet supercurrent 
flowing parallel to the interface of a simple F 1 /F 2 bilayer 
with noncollinear exchange directions, can significantly con¬ 
vert to a long-ranged spin-triplet supercurrent regardless of 
actual geometry. 

To further explore the generality of the phenomenon and 
provide an experimentally accessible platform to isolate it. 
Fig. 4 exhibits the diamagnetic supercurrent in a S/F1/F2/N 
structure [Fig. 1(c)] with c?_F2(=l-9Cs)^ > rffi(=0.1Cs) and 
<fjv=2.5£s, at three magnetization orientations of Fj: 
ai=n/2, /3i=0, 0.37T, and 0.57T, while /12 points along the 2 
direction; /32=0. As seen, Jy V flowing along the y direction, 
deeply penetrates F 2 laterally. Whereas, the singlet J*° and 
triplet Jy Z components are short-ranged and drop drastically 
midway through the magnet region. Hence, here is a clear and 
practical opportunity to generate extensive long-range spin- 
triplet supercurrents by parallel supercurrent flow relative to 
the interface of a uniformly magnetized Fj / F> structure with 
inequivalent F strips. The S/F1/F2/N structures proposed 
here are not only relatively simple to fabricate and readily 
accessible to experimental measurements [21], but they also 
serve as an efficient means to create samples for probing the 
previously discussed generic situations. In effect, the long- 
range triplet components deeply penetrate laterally into the 
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